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ABSTRACT 

We  continue  the  investigation  of  overcoming  Gibbs  phenomenon,  i.e.,  obtaining  exponen¬ 
tial  accuracy  at  all  points  including  at  the  discontinuities  themselves,  from  the  knowledge  of 
a  spectral  partial  sum  of  a  discontinuous  but  piecewise  analytic  function. 

We  show  that  if  we  are  given  the  first  N  expansion  coefficients  of  an  L%  function  f(x)  in 
terms  of  either  the  trigonometrical  polynomials  or  the  Chebyshev  or  Legendre  polynomials, 
we  can  construct  an  exponentially  convergent  approximation  to  the  point  values  of 
f(x)  in  any  sub-interval  in  which  it  is  analytic. 
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1  Introduction 


In  this  paper  we  continue  our  investigation  of  overcoming  Gibbs  phenomenon,  i.e.,  recovering 
pointwise  exponential  accuracy  at  all  points  including  at  the  discontinuities  themselves,  from 
the  knowledge  of  a  spectral  partial  sum  of  a  discontinuous  but  piecewise  analytic  function, 
which  we  started  in  [5]  and  [6]. 

Spectral  approximations  based  upon  trigonometric  polynomials  (Fourier,  for  periodic 
problems)  or  polynomials  (Chebyshev  or  Legendre,  for  non-periodic  problems)  are  expo¬ 
nentially  accurate  for  analytic  functions  [4],  [3].  However,  for  discontinuous  but  piecewise 
analytic  functions,  the  spectral  partial  sum  approximates  the  function  poorly  throughout 
the  domain.  Away  from  the  discontinuity  only  first  order  accuracy  is  achieved.  Near  the 
discontinuity  there  are  0(1)  oscillations  which  do  not  decrease  with  N ,  the  number  of  terms 
retained  in  the  spectral  sum.  This  is  known  as  Gibbs  phenomenon. 

In  [5]  and  [6]  we  treated  a  representative  problem  : 

•  Let  f(x)  be  an  analytic  but  nonperiodic  function.  Suppose  that  we  are  given  its  first 
—N  <k<  N  Fourier  coefficients.  Can  one  obtain  exponentially  convergent  pointwise 
approximations  in  the  maximum  norm? 

To  solve  this  problem  we  used  the  Gegenbauer  polynomials  CA(x),  which  are  orthogonal 
in  [—1, 1]  with  the  weight  function  (1  —  x3)A“&. 

The  procedure  of  overcoming  the  Gibbs  phenomenon  consists  of  two  steps: 

1.  Given  the  Fourier  partial  sum  of  the  first  N  terms,  we  first  recover  the  first 
m  ~  N  Gegenbauer  expansion  coefficients  with  exponential  accuracy.  This  can 
be  achieved  for  any  function,  as  long  as  we  choose  A  in  the  weight  function 
of  Gegenbauer  polynomials  to  be  proportional  to  N.  The  error  incurred  at  this 
stage  is  called  the  truncation  error. 

2.  The  next  step  is  to  prove,  for  an  analytic  function,  the  exponential  convergence 
of  its  Gegenbauer  expansion,  when  the  parameter  A  in  the  weight  function  is 
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proportional  to  the  number  of  terms  retained  in  the  expansion.  The  error  at  this 
stage  is  labelled  the  regularization  error. 

Thus  we  have  shown  how  to  overcome  the  Gibbs  phenomenon  for  an  analytic  but  non¬ 
periodic  function.  By  a  simple  shift  this  procedure  covers  the  case  of  any  analytic  function 
with  one  discontinuity. 

In  this  paper  we  treat  a  broader  class  of  problems.  Here  we  assume  that  f(x)  is  an  L2 
function  on  [—1, 1]  and  analytic  in  a  subinterval  [a,  6]  C  [—1, 1].  We  assume  that  the  spectral 
partial  sum  (based  on  either  the  Fourier  or  Chebyshev  or  Legendre  expansion)  of  a  function 
over  the  full  interval  [—1, 1]  is  known,  and  try  to  recover  exponentially  accurate  point  values 
over  a  subinterval  [a,  6]. 

We  will  follow  the  same  path  as  in  [5].  Basically  we  will  show  that  the  first  N  Fourier 
(or  Chebyshev  or  Legendre)  expansion  coefficients  contain  enough  information  such  that  a 
rapidly  converging  Gegenbauer  expansion  in  the  subinterval  [a,  6]  can  be  constructed.  As 
before,  we  will  separate  the  analysis  of  the  error  into  two  parts:  truncation  error  and  reg¬ 
ularization  error.  Truncation  error  measures  the  difference  between  the  exact  Gegenbauer 
coefficients  and  those  obtained  by  using  the  spectral  partial  sum.  These  will  be  investigated 
in  Section  3  for  the  Fourier  case  and  in  Section  4  for  the  Chebyshev  and  Legendre  cases. 
Regularization  error  then  measures  the  difference  between  the  Gegenbauer  expansion  using 
the  first  few  Gegenbauer  coefficients  and  the  function  itself  in  a  sub-interval  [a,  6],  in  which 
the  function  is  assumed  analytic.  The  estimation  of  the  regularization  error  is  in  Section  5. 
In  Section  2  we  will  provide  some  useful  properties  of  Gegenbauer  polynomials  to  be  used 
later.  Section  6  contains  a  summary  theorem  and  some  remarks. 

Throughout  this  paper,  we  will  use  A  to  denote  a  generic  constant  or  at  most  a  polynomial 
in  the  growing  parameters.  It  may  not  be  the  same  at  different  locations. 

The  message  of  this  paper  is  that  the  knowledge  of  the  partial  spectral  sum  of  an  L2 
function  in  [—1, 1]  furnishes  enough  information  so  that  an  exponentially  convergent  approx¬ 
imation  can  be  constructed  in  any  subinterval  in  which  f(x)  is  analytic. 
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In  particular,  if  /(«)  it  piecewiae  analytic,  then  an  exponentially  convergent  approxima¬ 
tion  in  the  maximum  norm  can  be  recovered  from  its  partial  spectral  expansion. 

2  Preliminaries 
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o 

We  will  seed  to  use  heavily  the  asymptotics  of  the  Gegenbauer  polynomials  for  large  n 
and  A.  For  this  we  need 
Lemma  2-S  (Stirling) 

For  any  number  x  such  that  i  >  1  we  have 

(2 *)***+*e-*  <  T(*  +  1)  <  (2ir)ix'+ie-eifc  (2.7) 

□ 


Lfmmi  3.4 

There  exists  a  constant  A  independent  of  A  and  n  such  that 

The  proof  follows  from  (2.6)  and  Stirling’s  formula  (2.7). 

□ 

In  the  analysis  we  will  extensively  use  the  relationship  between  the  Fourier  functions 
e****  and  the  Gegenbauer  polynomials  C*(x).  In  particular,  we  will  focus  our  attention  on 


the  Legendre  polynomials 

P.(z)  .  ct(z) 

(2.9) 

and  the  Chebyshev  polynomials 

T.(x)  =  jCJ(*) 

(2.10) 

Our  basic  formula  is  taken  from  [2],  page  213: 

Ltmmu  2iE 


^  />  -  ^-^imC,\x)dz  »  r(A)  (^)\'(1  +  A)J1+»(xn)  (2.11) 

where  J¥(x)  is  the  Bessel  function. 
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□ 


Lemma  2.5  gives  us  the  Gegenbauer  expansion  coefficients  of  the  Fourier  functions.  How¬ 
ever,  what  we  really  want  are  the  Fourier  coefficients  of  the  Gegenbauer  polynomials.  We 
have  those  coefficients,  luckily,  for  the  Legendre  and  Chebyshev  polynomials. 

kmmft  2,9 


Let  a?  be  the  Fourier  coefficients  of  the  Legendre  polynomial  Pn(x)  defined  in  (2.9)  and 
b%  be  the  Fourier  coefficients  of  the  Chebyshev  polynomial  Ts(x)  defined  in  (2.10),  i.e. 


ikirx 


and 


Pn(x)=  a%c' 

Jk=— oo 


r*(*)=  £  i?e“” 

fas- oo 


Then 


and 


*  m  IB1”*™ 


(2.12) 

(2.13) 

(2.14) 

(2.15) 

□ 


Equation  (2.14)  follows  from  (2.11)  upon  substituting  A  =  |  and  using  (2.6)  and  (2.3). 
Equation  (2.15)  can  be  found  in  [7],  page  836. 


□ 


Lemma  2.6  provides  us  with  a  better  understanding  of  the  terms  “small  scales”  and 
“large  scales”.  The  question  is  the  following:  it  is  clear  that  eiN*x  is  small  scale  for  large  N, 
intuitively  Ts(x)  and  Ps(x)  are  also  small  scale  functions,  however  the  expansions  (2.12) 
and  (2.13)  contain  all  the  scales.  So,  in  what  sense  are  Ts(x)  and  Pn(x)  small  scales? 

The  answer  is  given  in  the  following  Lemma: 

Lemma  2.7 


5 


Let  and  b£  be  the  Fourier  coefficients  of  the  Legendre  polynomial  Pjy(z)  and  the 
Chebyshev  polynomial  7jv(x)  respectively,  given  in  (2.12)-(2.15).  Then 


|aJ'|<Xmin 

lij'l  <  Amin  (l.tfjf)") 


(2.16) 


(2.17) 


where  A  is  a  constant  independent  of  k  or  N. 


ProQfc 


We  start  with  the  following  estimate  for  the  Bessel  function  Jn{nz)  (see  [1],  page  362): 


|J„(nz)| 


A  Y 


0<z<  1 


(2.18) 


From  (2.18)  we  can  deduce  trivially  that 


M.MI  <  ( j)"  o<2<i 


We  also  have  the  uniform  bound  on  Bessel  functions 


(2.19) 


\U*)\  <  1 


(2.20) 


Using  (2.19)  and  (2.20)  in  (2.14)  and  (2.15),  with  z  =  $  and  n  =  N,  we  obtain  (2.16) 
and  (2.17). 


Lemma  2.7  provides  the  answer  why  Pn{x)  and  Tjv(x)  (for  large  N)  can  be  considered 
as  small  scale  functions.  Basically  only  the  terms  that  are  greater  than  N  appear  in  their 
Fourier  expansions  (2.12)  and  (2.13).  The  lower  terms  decay  exponentially  with  N. 
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3  Truncation  Error  in  a  Sub-interval — Fourier  Method 


Consider  an  arbitrary  Lj  function  f(x)  defined  in  [—1,1].  Suppose  that  the  first  2N  +  1 
Fourier  coefficients  of  this  function  are  given 


m  =  i  £  i*i  <  n 


In  [5]  we  showed  that  we  can  recover  (with  exponential  accuracy)  the  first  m  ~  N  terms 
in  the  Gegenbauer  expansion  of  /(x)  based  on  the  interval  [—1,1]. 

Here  we  are  interested  to  find  the  Gegenbauer  expansion  of  /(x)  based  on  a  sub-interval 
[o,  b]  C  [—1,1].  This  will  allow  us  to  handle  multiple  discontinuities  and  discontinuities  of 
unknown  locations.  We  will  show  that  in  this  case  too,  we  can  get  exponential  accuracy.  We 
start  by  introducing  the  local  variable  £: 


The  local  variable  £  is  defined  by 


x  =  x(£)  =  e£  +  6 


where 


b—  a 


b  +  a 


Thus  when  a  <  x  <  b,  —  1  <  £  <  1. 


Denote  now  the  Fourier  partial  sum  by  //v(x),  namely 

fs(x)  =  £  f(k)eik*x  (3.4) 

k=-N 

As  mentioned  before,  we  are  trying  to  recover  the  Gegenbauer  expansion  coefficients 
based  on  the  sub-interval  [o,6],  i.e.  we  would  like  to  find  the  first  m  ~  N  coefficients  /(A(/) 


in  the  expansion 


/(*«))  =  EKwH  0 

1=0 
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where 


mi)  =  /_'(!  -  (3.6) 

Of  course  we  do  not  have  /eA(/)  at  our  disposal  but  only  an  approximation  based  on  the 
Fourier  partial  sum  //v(x),  thus  we  have 

j,A(<)  =  £(  1  -  (3.7) 

The  Truncation  Error  is  the  measure  of  deviation  of  the  approximate  Gegenbauer  expan¬ 
sion  from  the  true  expansion: 

TE(  A,  m,  JV, «)  =  max  |  f(/‘(/)  -  jtA(i))C,A({)l  (3.8) 

/=0 

We  can  show  that  the  truncation  error  is  exponentially  small: 

Theorem  3,3; 

Let  /(x)  be  an  L2  function  whose  first  27V -fl  Fourier  coefficients  f(k)  defined  in  (3.1)  are 
known.  Let  g*(l)  be  the  Gegenbauer  coefficients,  based  on  the  sub-interval  [a,  6]  expansion 

A 

of  the  Fourier  partial  sum  /n{x)  defined  in  (3.7),  and  ft(l)  be  those  of  /(x)  defined  in  (3.6). 
Assume  that 

A  =  aeN  m  =  0eN  (3.9) 

Where  e  is  the  half  length  of  the  sub-interval  [a,  6]  defined  in  (3.3).  Then  the  truncation 
error  is  bounded  by 

TE{aeN ,  0tN,  N,  e)  <  AqtN ,  (3.10) 


where 


_  {0  +  2 a)fi+2a 

q~  (2*e)a<Xa0fi 

and  A  grows  at  most  quadratically  with  N.  In  particular,  if  a  =  0  «  A,  then 


(3.11) 


q  «  0.8 


Prevf; 
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The  proof  follows  exactly  the  one  in  [5],  Theorem  3.1  and  3.3,  with  proper  accounting 
for  the  interval  length  c. 

□ 

4  Truncation  Error  in  a  Sub-interval — Chebyshev  and 
Legendre  Methods 

In  this  section  we  assume  that  we  have  the  first  N  +  1  Chebyshev  or  Legendre  coefficients 
in  the  expansion  of  an  La  function  f(x)  based  on  the  interval  —  1  <  x  <  1.  We  show  how 
to  get  its  Gegenbauer  expansion  based  on  a  sub-interval  [a,  b]  C  [—1,1]  with  exponential 
accuracy  in  N.  The  proof  makes  use  of  Lemma  2.7,  which  states  that  the  low  modes  of  the 
Fourier  expansion  of  the  Chebyshev  or  Legendre  polynomial  of  degree  N  decay  exponentially 
with  N.  Thus  only  the  high  modes  count,  and  those  are  shown  in  Section  2  to  have  small 
Gegenbauer  coefficients. 

For  simplicity  of  notation  we  will  deal  with  c\  (x),  which  is  the  Legendre  polynomial 
P*(x)  (see  (2.9)),  and  Cj?(x),  which  is  £  times  the  Chebyshev  polynomial  T*(x)  (see  (2.10)). 
The  Chebyshev  or  Legendre  expansion  coefficients  of  /(x)  are  given  by 

/'(*)  =  ^  jj  i  -  x'y-if(x)c;(z)dz  (4.i) 

for  n  =  0  or  p  =  |. 

In  this  section  we  consider  functions  /(x)  satisfying 
Assumption  4.1: 

|/'*(Jk)|  <  A  for  all  ifc. 

□ 

We  remark  that  if  /(x)  is  an  Li  function  the  assumption  is  fulfilled. 

We  assume  that  we  know  the  first  N  +  1  Chebyshev  or  Legendre  coefficients,  /**(&)  for 
0  <  k  <  N,  and  define 

f^x)  =  'tf“(k)C;(Z)  (4.2) 
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Note  that  /£(x)  does  not  converge  fast  to  f(x)  if  there  exist  discontinuities  inside  the  domain. 

The  function  f(x)  has  also  a  Gegenbauer  expansion  in  a  sub-interval  [a,  6].  With  £,  e  and 
6  defined  in  (3.2)-(3.3),  we  have 

m  +  S)  =  tf\‘)C,y(0  (4.3) 

IssO 

where  the  Gegenbauer  coefficients  /*(/)  are  defined  by 

/«'(')  =  ^  £(1  -  nA-i/(<{  +  (4.4) 

As  before,  we  do  not  have  /*(/)  at  our  disposal,  but  only  an  approximation  based  on  the 
partial  Chebyshev  or  Legendre  sum  /Jv(x),  thus  we  have 

ifw = A  jji  -  + i)cHm  (4.5) 

How  well  do  g*(l)  approximate  /A(/)?  To  answer  this  question  we  define 
Definition  4.2: 

The  truncation  error  is  defined  by 

TE(X,  m,  1 V,  e)  =  max  |  f  )(/»(0  -  j,A(OK?(OI  (4.6) 

Is 0 

where  /A(/)  are  defined  by  (4.4)  and  <jA(/)  are  defined  by  (4.5). 

O 

The  truncation  error  is  the  measure  of  the  distance  between  the  true  Gegenbauer  expan¬ 
sion  in  the  interval  [a,  6]  and  its  approximation  based  on  the  Chebyshev  or  Legendre  partial 
sum  in  [—1, 1]. 

We  first  have  the  following  lemma: 

Lemma  4.3: 


TE(\,m,N,e)<  £  |/“(,)|£ 

g=N+ 1  1=0 


/->(1 "  W  +  «)c,A(f)<if 


(4.7) 
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Proof: 

From  (4.4)  and  (4.5)  we  have 

|/A(0  -  «.A(i)|  =  jjr  |/_V  -  «*)*-»(/(«{  +  «)- /»«  +  «))C(A«)d{ 

Substituting  (4.8)  into  (4.6),  recalling  (2.4)  and 

/(«*+*)-#(*+*)«  £  rwfct+s) 

q=N+l 

we  obtain  (4.7). 


(4.8) 

(4.9) 

□ 


For  simplicity  of  notations  we  denote: 


^r,U  = 


hr 


9>  JV 


(4.10) 


In  order  to  estimate  this  term,  we  consider  the  Fourier  series  approximation  of  C*(e{  +  6) 
from  (2.12)  and  (2.13)  to  get 

4  4: 


F,.U  <  ^C,J(l)r(A)  g  (j^)A  |J,+A(farc)|  |«J|  (4.11) 

where  a9k  is  given  either  by  (2.14)  or  by  (2.15),  and  A  grows  at  most  as  fast  as  a  linear 
polynomial  in  A  and  /. 

Proof: 

We  expand  C*(x)  using  (2.12)  or  (2.13): 


OO 


Cq{x)  =  Y  oj^e’***  =  Y 

fc=— oo  k=—oo 

where  a*  ^  =  a*  is  given  by  (2.14)  and  a*0  =  is  given  by  (2.15). 

We  now  use  (2.11)  to  get 

=  rw(jj^)  (f+Wwt*™)! 


(4.12) 


(4.13) 


Substituting  (4.12)  and  (4.13)  into  (4.10)  gives  us  (4.11). 
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□ 


We  are  ready  now  to  estimate  directly  F9jt\  defined  in  (4.10). 

Lemma  4.5: 

Consider  defined  in  (4.10).  Let  A  =  ^€JV  and  l  =  ~f\  where  7  <  1.  Then  for  every 
q>  N  we  have 


F,,i,a  <  Ae~x 


(?)  - 


(4.14) 


where  A  grows  at  most  quadratically  with  N. 

Proof: 

Let  kq  —  be  the  integer  part  of  We  split  the  sum  on  the  right  hand  side  of  (4.11) 
into 

f,m<a(z+  E  )  b(i)r(A)  Ml+»(fcr«)l  Ml]  =  A  (ff&  +  F%)  (4.15) 

\Asl  L 

We  first  use  (2.3)  and  (2.7)  to  obtain: 

«,IS.(«I±2=)‘S.(2)‘ 


(4.16) 


and 


(4.17) 


T(A)  <  AAV* 

We  start  with  an  estimation  of  the  second  sum 

*=fcf+ 1 

We  use  (4.16),  (4.17),  the  fact  that  |«7j+a(£)1  <  1  and  |o||  <  1,  and  k  >  kq  =  [^]  and 
A  =  to  get 


flZ=  E  c>(i)r(x)(£)xVMii« 


^-(7)  ■ 


(4-18) 


where  A  grows  at  most  linearly  with  N,  the  growth  coming  from  the  summation.  For  the 
first  sum 

=  E  c*(i)r(A)  (£)  |a+a(*t£)|  mi  , 
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we  still  use  (4.16),  (4.17)  and  the  bound  |Ji+;t(x)|  <  1.  However,  for  the  estimate  on  |aj[|, 
we  use  (2.16),  (2.17): 


V2 


and  the  fact  q  >  N  >  A  to  get 


and,  upon  substituting  for  kq  =  and  A  =  j|_c7V: 

Combining  (4.18)  and  (4.21)  we  obtain  (4.14). 


(4.19) 


(4.20) 


(4.21) 


We  are  now  ready  for  the  main  theorem  of  this  section: 

Throrani  4t6; 

Let  the  truncation  error  be  defined  in  (4.6).  Let  A  =  ^-N  and  m  <  A,  then 

TE{ A,  m,  JV,  e)  <  i4e_A  (4.22) 

where  A  grows  at  most  as  a  fourth  order  polynomial  in  N . 

Proof; 

The  theorem  follows  from  (4.7),  the  Assumption  4.1,  and  (4.14). 

□ 


5  Regularization  Error  for  a  Sub-interval 

In  this  section  we  study  the  second  part  of  the  error,  which  is  caused  by  using  a  finite 
Gegenbauer  expansion  based  on  a  sub-interval  [a,  6]  C  [—1,1],  to  approximate  a  function 
f(x)  which  is  assumed  analytic  in  this  sub-interval  [a,  b ].  Since  in  previous  sections  we  have 
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established  exponentially  small  truncation  errors  if  A  and  m  are  both  growing  linearly  with 
N,  in  this  section  we  will  consider  the  case  of  A  ~  m. 

We  will  assume  that  f(x)  is  an  analytic  function  cm  [a,  b]  satisfying 

Assumption  5.1 

There  exists  constants  p  >  1  and  C{p)  such  that,  for  every  ifc  >  0, 

<  (s.i) 

□ 


max 

a<*<fr 


This  is  a  standard  assumption  for  analytic  functions,  p  is  the  distance  from  [a,  6]  to  the 
nearest  singularity  of  f(x)  in  the  complex  plane  (see,  for  example,  [8]).  This  assumption  and 
the  following  proof  can  be  modified  to  allow  for  an  arbitrary  ellipse  region  around  [a,  b]  for 
the  analytidty  of  f(x)  (H.  Vandeven,  private  communication). 

Let  us  consider  the  Gegenbauer  partial  sum  of  the  first  m  terms  for  the  function  /(e£+6): 


/£*«)= £/,*(')(?(() 

IsO 


(5.2) 


with  (,  e  and  6  defined  by  (3.2)  and  (3.3),  and  the  Gegenbauer  coefficients  based  on  [a,b] 
defined  by 

/‘(i)  =  A  Jj  1  -  e’H/K + WJ«)<f£  (5.3) 

We  want  to  estimate  the  regularization  error  in  the  maximum  norm: 


RE(X,m,e)=  max 


/K+«)-E£(')c?(() 

1=0 


(5.4) 


We  have  the  following  result  for  the  estimation  of  the  regularization  error,  when  A  ~  m: 

Theorem 


Assume  A  =  7 m  where  7  is  a  positive  constant.  If  f(x)  is  analytic  in  [0,6]  C  [—1,1] 
satisfying  the  Assumption  5.1,  then  the  regularization  error  defined  in  (5.4)  can  be  bounded 
by 

REfam,  m,  e)  <  Aqm  (5.5) 
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where  q  is  given  by 


9  p21+^Y»(l+7)1'H  1  ; 

which  is  always  less  than  1.  In  particular,  if  7  =  1  and  m  —  (3N  where  fl  is  a  positive 


constant,  then 


RE(/3N,  0N,  e)  <  AqN 


=  (-V 


Ptftgfc 


The  proof  follows  in  exactly  in  the  same  fashion  as  in  [5],  Lemma  4.2,  Theorems  4.3  and 
4.4,  with  proper  accounting  for  the  transformation  x(£)  =  +  S. 


6  Concluding  Remarks 

In  this  Section,  we  combine  the  estimates  for  truncation  errors  and  regularization  errors  in 
previous  sections,  to  obtain  the  main  theorem  of  this  paper: 


e  Gibbs 


Consider  an  L%  function  /(x)  on  [—1, 1],  which  is  analytic  in  a  sub-interval  [a,  b]  C  [—1, 1] 


and  satiates 


'-1 


Assume  that  either  the  first  2 N  +  1  Fourier  coefficients 


/(*)  =  X-  Jl  f(x)e~ikir*dx ,  -N  <  k<  N 


or  the  first  N  +  1  Chebyshev  or  Legendre  coefficients 


/'(*)  =  jj  1  -  xJf-»/(x)Cr(*)<fa,  „  =  0  or  „  =  i 
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are  known.  Let  g,(l),  0  <  /  <  m  be  the  Gegenbauer  expansion  coefficients,  based  on  the  sub¬ 
interval  [a,  b],  of  the  spectral  partial  sum  fs (*)  in  (3.4)  or  /£(*)  in  (4.2).  These  Gegenbauer 
coefficients  are  defined  in  (3.7)  for  the  Fourier  case  and  in  (4.5)  for  the  Chebyshev  or  Legendre 
case.  Thai  for  A  =  m  =  fteN  where  ft  <  (Fourier  case)  or  ft  <  (Chebyshev  or 
Legendre  case),  we  have 

m +«)- <  a (4w + &)  (6.2) 

lm 0 


max 


where 


for  the  Fourier  case  and 


for  the  Chebyshev  or  Legendre  case,  and  A  is  at  most  a  fourth  order  polynomial  in  N. 

Proof? 

Just  combine  the  results  of  Theorems  3.2,  4.6  and  5.2. 


□ 


We  finally  give  two  remarks: 

Remark  g«3 

It  can  be  seen  from  (6.2/  tuat  the  truncation  error  gets  bigger  when  the  interval  half- 
length  e  decreases,  because  the  parameters  m  and  A  must  decrease  with  t.  This  in  turn  also 
affects  the  regularization  error,  although  there  is  a  factor  e  in  which  offsets  this  effect. 

□ 


No  attempt  has  been  made  in  thi*.  to  optimize  the  parameters. 


□ 
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